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Abstract. In recent years, a number of extragalactic massive star clusters 
that are still deeply embedded in their birth material have been discovered. 
These objects represent the youngest stage of massive star cluster evolution 
yet observed, and the most massive and dense of these may be proto globular 
clusters. Their properties appear to be similar to those of ultracompact HII 
regions in the Galaxy, but scaled up in total mass and luminosity. In many 
cases, these clusters are only visible at mid-IR to radio wavelengths, and they 
have typically been detected as "inverted" spectrum radio sources. However, the 
set of existing observations is anemic, and our current physical model for these 
natal clusters in simplistic. This article will overview what we think we know 
about these objects based on existing observations and outline some of the most 
significant gaps in our current understanding. 



1. Introduction 

There is no question that massive star clusters, and super star clusters in partic- 
ular, are among the most extreme star formation environments in the local uni- 
verse. Moreover, these clusters can have a dramatic affect on their surrounding 
interstellar medium, and in some cases the intergalactic medium. Nevertheless, 
the conditions required for their formation remain far from understood. As just 
one example, there is still significant debate over whether super star clusters 
(SSCs) are merely the result of populating a cluster initial mass function statis- 
tically, or whether they truly require a separate mode of star formation. While a 
multitude of adolescent massive star clusters have been observed at optical wave- 
lengths, these observations cannot penetrate the natal material surrounding the 
youngest clusters in order to reveal the properties of their birth environments 
(Figure^). Nevertheless, if we want to understand the formation of massive star 
clusters, it is not a terribly bad idea to observe them while they are forming. 

A number of investigations at infrared (IR) and radio wavelengths have 
now revealed a populati on of extremely young; ma s sive star clusters in a growing 
sample of galaxies fe.g. Ilurner. Ho. Sz Becklll998l: iKobulnickv fc Johnsonlll99 " 



Tarc hi et ali 2000; Beck. Turner, fc Gorii an 2001; Goriian. T urner, fc Beck 20 0l|; 
Vacca. Johnson, fc Gontil2002l:I.Tohnson et al.l200ll:lNeff fc Ulvestadl2000l:lPlante fc Sauvagel 



2002| : lBeck et al.l2002l : lJohnson. Indebetouw. fc Pisanol2003l : I Johnson fc Kobulnickvl " 
2003). These observations suggest that the early stages of massive star clus- 
ter evolution may parallel those of individual massive stars observed in the 
Milky Way. Therefore, our knowledge about individual ultracompact Hll regions 
(UCHns) is critical to understanding and interpreting observations of embedded 



1 



2 



Johnson 



Evolutionary State 




Observable Regimes 


Deeply Embedded 
Ultra-young Cluster 




Mid-IR - Radio (cm) 


Highly Extincted 
Young Cluster 


T 


Near-IR- Radio (cm) 


Fully Emerged 
Adolescent Cluster 


- ' ; ' : 


Optical - Near-IR 



Figure 1. A cartoon illustrating stages in the early evolution of massive 
star clusters. Observations in at thermal infrared to radio wavelengths are 
required to probe their birth environments. 



massive star clusters. Of course, individual massive star formation is also poorly 
understood, and we would be remiss to blindly extrapolate the inferred prop- 
erties of individual massive proto-stars to the formation of SSCs. Perhaps the 
most obvious way in which this extrapolation could lead us astray is by neglect- 
ing the effects of the density of proto-stars on star formation. In fact, it seems 
likely that a tremendous amount of physical insight will come from trying to 
understand how the properties of star formation scale between the two extremes 
of individual massive star formation and SSC formation. 

Despite the pioneering observations done by a number of investigators, our 
current knowledge about embedded massive star clusters is anemic, and we have 
only begun to investigate their properties, environments, and evolution. In the 
pages that follow, I will discuss the current observational limitations (0, give 
a brief tour of the sample of extragalactic embedded clusters (3&Q), overview 
what we think we know about these objects (as well as what we do not know, 
tj4.|) . and outline the impact that the next generation of powerful observatories 
at infrared to radio wavelengths will have on this area of research fflo~|). 

2. Observational Limits 

Because star formation takes place within dense and dusty cocoons, we can 
only probe these regions with observations at infrared or longer wavelengths. In 
fact, some natal clusters (presumably the you ngest) are even invisible at near -IR 
wavelengths as long as the K-band (~ 2.2^m. lVacca. Johnson, Contill2002l ) or 
even the L-band (~ 3.8/im, M. Sauvage, private communication). On the other 
hand, radio emission is completely unfettered by dust, but the observations need 
to be sensitive at high frequencies (> 5 GHz) in order to avoid self-absorption by 
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Figure 2. The 5<r sensitivities for 8 hr integrations at 6 cm for various radio 
interferometers (both current and planned) are indicated by the horizontal 
grey shaded areas. The flux densities of "typical" UCHIIs and natal super 
star clusters as a function of distance are shown with the diagonal lines. The 
distances to a handful of reference galaxies are indicated. 



the dense gas and disentangle non-thermal radio emission. Of course, obtaining 
the necessary sensitivity is only half of the battle; achieving adequate spatial 
resolution is also a limiting constraint with current facilities. Roughly speaking, 
a linear resolution of at least ~ 0.5 pc is necessary to observe individual UCHns, 
and embedded SSCs require a linear resolution of ~ 10 pc or better. 

In principle, mid-IR to sub-millimeter wavelengths are ideal for observing 
extragalactic UCHns and natal SSCs because these objects are extremely lumi- 
nous in this regime due to thermal dust emission (Figure [SJ. As a demonstration 
of their potency at thermal wavelengths, UCHns are among the most luminous 
Galactic objects in the IRAS catalog ( > 60% of the IRA S sources with flux den- 
sities > 10 4 Jy at lOOjum are UCHiis. IWood fc Churchwel]lll98<l . and the em- 
bedded SSCs in the starburst galaxies He 2-10 and NG C 5253 are responsible for 
nearly all of the mid-IR emission from these galaxies dBeck. Turner. Goriianl 
200l|; iGoriian. Turner, fc Beckll200ll : IVacca. Johnson, fc Contill200^ . However, 
we are severely limited by spatial resolution with the facilities currently available 
at infrared to sub-millimeter wavelengths. This situation will improve to some 
degree with the next generation of observatories that are beginning to come on- 
line; however, even the Spitzer, Herschel, and SOFIA telescopes will not obtain 
spatial resolutions of better than ~ 1" at their shortest wavelengths. 

The sensitivity and spatial resolution of existing radio interferometers is 
somewhat better matched to observing extragalactic UCHiis and natal SSCs. 
Figure |2] illustrates the sensitivities of various radio interferometers. The cur- 
rent facilities have the sensitivity to detect individual UCHiis nearly as far away 
as the starburst galaxy NGC 253, or embedded SSCs nearly out to ultralu- 
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Figure 3. Examples of the most nearby and most distant embedded clus- 
ters yet studied, (left) The Hn region N81 in the SMC; an HST Ha im- 
age is shown in gray-scale, and ATCA 3 cm radio emission is shown in con- 
tour (In debetouw. Johnson, fc Contil 120041) . (right) The starburst galaxy 
SBS0335-052; an ffSTI-band image is shown in gray-scale, and VLA 1.3 cm 
radio emission is shown in contour (Johnson & Plante in prep.) 



minous infrared galaxy Arp 220. Perhaps the most obvious places to look for 
extragalactic UCHns are the Magellanic Clouds, however the best spatial resolu- 
tion available in the southern hemisphere from the Australia Telescope Compact 
Array at this distance is ~ 1 pc, which limits our ability to asses whether Hn re- 
gions are actually "ultra" compact. Another obvious place to search of UCHns 
is in M33, the sister galaxy to the Milky Way. However, a linear resolution 
of better than ~ 0.1" is required in order to observe UCHns at this distance 
(840 kpc), which will be difficult (but not impossible) to achieve with the Very 
Large Array (VLA). The constraints are somewhat more relaxed for observing 
embedded SSCs, which are both significantly larger and brighter than individual 
UCHiis. With the VLA, it is possible to resolve embedded SSCs in galaxies out 
to distances of ~ 10 Mpc. Beyond this distance, both crowding and contamina- 
tion by non-thermal and optically-thin thermal radio emission becomes an issue, 
and one has to use caution when interpreting the data. 



3. Sample of Extragalactic Natal Clusters 

A number of embedded massive star clusters have been identified in a variety 
of galaxies. The objects that have been most well-studied to date are those 
that were identified in the incipient radio observati ons of the starburst galax- 
ies NGC 5253 (~ 4 Mpc) an d He 2-10 (~ 9 Mpc) (|Turner. Ho. fc Becldll99Sl : 
iKobulnickv Johnsonl ll999). Subsequent work has revealed a number of em- 
bedded star form ing regions in a range of galact i c env ironments. In the most 
nearby examples, llndebetouw. Johnson, h Contil ((20041 ) have identified a num- 
ber of candidate UCHiis and UCHu complexes in the Magellanic Clouds (here 
the word "candidate" is used because the spatial resolution is not adequate to 
determine whether the regions are "ultra" compact). In this sample, the radio 
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observations indicate stellar ionizing fluxes and morphologies consistent with 
sources ranging from single UCHiis to relatively modest UCHu associations 
(Figure Oh)- 

Compact radio-detected clusters have also been identified in galaxies such 
as NGC 6822 at ~ 0.5 Mpc, IC 4662 at ~ 2 Mpc, IIZw40 at ~ 11 Mpc, 
Haro 3 at ~ 13 Mpc, NGC 5398 at ~ 17 Mpc, NGC 214 6 at ~ 15 Mpc 
the Antennae at ~ 21 Mpc, and SBS 0335-052 at ~ 53 Mpc dNeff fc Ulyestadl 
20001: iTarchi et al.ll2000l : iBeck et alJl2002l : I.Tohnson. Indebetouw. fc Pisanoll2003T 
Johnson et al.l 120041 . Johnson & Plante in prep., Seth & Johnson in prep.) . 
The star forming regions in this sample have properties ranging from those of 
natal OB-associations to incredibly massive SSCs. For example, the main radio 
source in SBS 0335-052 has an ionizing flux equivalent to ~ 12, 000 07-type 
stars (Figure OJa; note that the linear resolution of these observations is only 
~ 75 pc, and this source is likely to be composed of several individual clusters). 

4. What We Think We Know & What We Know We Don't Know 

Given the limited observations available at the present time, attempts to con- 
strain the physical properties of embedded clusters have been fairly crude. As a 
result, the current physical model we have for these objects is simplistic: When 
an SSC is born, the massive stars ionize the surrounding natal material, creat- 
ing an extremely dense Hn region. The dense Hn region is in turn surrounded 
by a dust cocoon that is presumably quite warm on the inner boundary, with 
the temperature decreasing at distances farther from the ionizing stars. The 
constituent ionizing stars may be surrounded by individual Hn regions and dust 
cocoons in the earliest stages, but this is currently unknown. Clusters that are 
more sparse are likely to resemble ultracompact Hn region (UCHu) complexes 
in the Milky Way (such as W49A), and contain spatially discrete Hn regions 
and cocoons, but the existence of individual cocoons becomes less tenable in the 
extreme stellar densities of the cores of massive SSCs. 

The dense Hn regions are characterized by an "inverted" radio spectral 
energy distributions due to optically thick free-free emission (a > 0, where 
S u oc z/ a ). The frequency at which the radio emission becomes optically thick 
is largely dependent on the density of the Hn region. In the earliest stages of 
cluster evolution, all of the stellar luminosity is reprocessed by the surrounding 
dust cocoon, and this cocoon is observable at thermal infrared to sub-millimeter 
wavelengths. To some extent, the nature of the resulting spectral energy distri- 
bution depends on whether the dust in the cocoon is smooth or clumpy, but this 
is not currently constrained by observations. 

Existing observations have provided estimates for some of the physical prop- 
erties of natal clusters, including size, electron density, pressure, stellar mass, 
Hn mass, dust mass, and age. For the most nearby clusters, their size can 
be directly measured using high spatial resolution observations. For more dis- 
tant clusters, model Hn regions can be used to fit the radio spectral energy 
distribution and infer both the radius and electron density (Figure HJd) . The 
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Figure 4. (left) Natal SSCs in the starburst galaxy He 2-10. An HST 
I-band image is shown in gray-scale, and VLA 2cm emission is shown in 
contour. The four radio sources without optical counterparts are natal SSCs. 
The arrow indicates the source for which the spectral energy distribution is 
shown, (right) The radio fluxes for the most luminous natal SSC in He 2-10 
along with the best fit spectral energy distribution for a model Hn region. 



embedded SSCs that have b een analyzed in this way (e.g. those in He 2-10, 
Uohnson fc Kobulnickvl 1200.1 ) have Hn region radii of only a few parsecs, and 
the electron densities have global values of ~ 10 3 cm" 3 , with peak values in 
excess of ~ 10 6 cm -3 . Assuming temperatures of ~ 10 4 K, these densities imply 
global pressures of P/ks > 10 7 K cm -3 and peak pressures reaching values of 
P/k B ~ 10 10 K cm" 3 . 

The masses associated with the stellar clusters, Hn regions, and dust co- 
coons can be bootstrapped from radio and infrared observations. The optically- 
thin thermal flux (typically measured with the highest available radio frequency) 
can be used to infer a total ionizing luminosity, and the ionizing luminosity 
can be translated into an embedded stellar mass by assuming a stellar initial 
mass function. Using this technique, clusters with a range of masses have been 
detected (spanning modest OB-associations to massive SSCs). The Hn mass 
can be estimated given the radii and densities for the ionized regions. In at 
least some cases, the Hll mass has anomalously low values (< 5% of the stellar 
mass, Uohnson fc Kobulnickvl 12003 ), which we tentatively interpret as a sign of 
the youth of the Hn regions. The mass of the dust cocoons around embed- 
ded SSCs has also been estimated in a few cases using infrared observations; 
these estimates have suggested dust masses of ~ 10 5 — 10 6 M^ for massive SSCs 
in SB S0335-052 and He 2-10 (|Plante fc Sauvagdl2002l : IVacca. Johnson, fc Contil 
I2002D . 

The ages of the natal clusters are exceptionally tricky to ascertain, as there 
are no available spectral diagnostics to determine age at wavelengths longer than 
the near-IR (and for shorter wavelengths the diagnostics are unreliable for ages 
less than ~ 1 Myr in any case). The best we can do is offer circumstantial 
evidence and look at the back of a proverbial envelope. There are two inde- 
pendent "back-of-the-envelope" calculations that can be used to constrain the 
lifetimes of the embedded phase of SSC evolution: (1) The pressures inferred for 
these regions are several orders of magnitude higher than those typical in the 
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ISM. As a result, the Hn regions ought to expand and dissipate on very short 
time-scales (r < 10 6 yrs). However, there are caveats to this argument: First, 
we do not know the pressure of the surrounding ISM, and it could be a good 
deal higher than normal since it is a starburst region. Second, if the clusters are 
massive enough, the gravity could be large eno ugh to help balance pressure and 
confine the Hn region (e.g. iTurner et alJ 120031 ). (2) If star formation has been 
relatively continuous over ~ 10 Myr, then a comparison between the total num- 
ber of clusters with ages < 10 Myr and the number embedded clusters ought to 
reflect the their relative lifetimes. W hen this has been applied to galaxies with 
good enough statistics (e.g. He 2-10. iKobulnickv h Johnsonlll999h . the clusters 
appear to spend ~ 10 — 20% of the lifetimes of the constituent massive stars 
in the embedded phase (r ~ 0.5 — 1 Myr). There are caveats to this argument 
as well. First, it relies on star formation being continuous over the lifetimes of 
massive stars. Second, we have to make sure that both the optical and radio 
observations are sensitive to roughly the same cluster masses. In any case, it is 
comforting that both calculations suggest the same lifetimes for the embedded 
phase of cluster evolution, so our estimates are probably not too far off. 

While we have estimates for a few of the quantities associated with the 
birth of massive star clusters, there are a potpourri of properties that it would 
be nice to know more about. For example, it would be instructive to directly 
measure densities, pressures, and temperatures using various line diagnostics at 
IR to radio wavelengths. Another quantity that might be particularly useful 
for theorists (that is almost completely unconstrained) is the star formation 
efficiency in these objects. We also have no idea how much ionizing radiation 
could be escaping from the natal cocoons (and it could be quite a lot if the 
dust is clumpy). In fact, we know almost nothing about the natal dust cocoons 
(including the temperatures, grain compositions, and geometry). This list is 
certainly not exhaustive. An understanding of these quantities, and how they 
relate to the properties of the resulting cluster, will eventually provide a great 
deal of insight into the requirements for massive star cluster formation. 



5. Future Prospects 

While there are no lack of questions remaining about the birth of massive star 
clusters, it is an exceptional good time to be asking these questions. Over the 
next decade (or so), a range of powerful facilities are going to become available 
at infrared to radio wavelengths. Figure |S] illustrates the spatial resolution and 
sensitivity of the IR to millimeter facilities that will soon become available. In 
the thermal infrared, Spitzer, Herschel, and SOFIA will provide a vast improve- 
ment over previous facilities, and embedded star forming regions are likely to 
be a dominate source of flux in the thermal spectral energy distribution of star 
forming galaxies. However, we will have to be cautious in our interpretation 
of data from these observatories because they will not have the spatial reso- 
lution necessary to isolate individual natal super star clusters. At millimeter 
wavelengths, the Atacama Large Millimeter Array will open up an entirely new 
window for observing extragalactic star formation with resolution and sensitiv- 
ity incredibly well- matched to observing embedded clusters. In the radio regime 
(see Figure 12), the Expanded Very Large Array and Square Kilometer Array are 
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Figure 5. (left) The anticipated spatial resolution and wavelength coverage 
for upcoming infrared and millimeter observatories, (right) The sensitivities 
of these observatories normalized to 500s integration time. A model spectral 
energy distribution (Whitney et al. in prep) for a 10 6 M© natal SSC at a 
distance of 10 Mpc is over-plotted with a dashed line for reference. 



going to achieve sensitivities and resolutions that will allow us to study nearby 
natal clusters in fantastic detail and enable us to discover new natal clusters out 
to distances of 100 Mpc and beyond. 
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